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Evaluation of three- and four-point diagrams with massless internal particles and 
arbitrary external momenta is considered. Exact results for some two-loop diagrams 
(planar and non-planar three-point contributions and the "double box" diagram) are 
obtained in terms of polylogarithms. 

1. Recent development of modern accelerators requires information about 
higher-order loop corrections to elementary particle processes. Some of the in- 
teresting problems (corrections to multijet processes, Bhabha scattering, etc.) are 
connected with the evaluation of two-loop three- and four-point diagrams. The 
case of massless internal propagators is of special interest, because often we are 
confronted either with really massless particles (photon, gluon) or with particles 
whose masses can be neglected in high-energy processes (electron, light quarks, 
etc.). For example, in a number of publications 0, ||, [| (see also [Q and references 
therein) the asymptotic behaviour of ladder diagrams has been examined for high 
energies and momentum transfer in leading logarithmic approximation. 

The present paper is devoted to some methods which enable one to obtain exact 
expressions for some types of Feynman diagrams with massless internal particles 
and arbitrary external momenta (with three and four external particles). The used 
approach involves the following tools: (i) the Feynman parametric representation, 
(ii) the "uniqueness" conditions (see, e.g., in H), (iii) Mellin-Barnes contour inte- 
grals and (iv) Fourier transform to the coordinate space. We shall consider only 
scalar diagrams (corresponding to massless <f) 3 theory), because expressions occur- 
ring in realistic calculations can be reduced to such scalar integrals (see, e.g., H). 

The remainder of the paper is organized as follows. In Section 2 we present the 
main steps of the approach via examples of one-loop triangle and box diagrams. 
In Section 3 we apply this technique to obtain expressions for two-loop three- and 
four-point ladder (planar) diagrams. Section 4 discusses the evaluation of some 
non-planar graphs. Section 5 formulates and discusses the main results. 

2. Let us start by some definitions and results for massless triangle diagrams 
(see Fig. 1). Here and below we shall consider all external momenta to be ingoing 
(.Pi +P2 +P3 = 0). The corresponding Feynman integral is 
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Figure 1: 



where q 3 — q-x = pi, qi — q 3 = P2, q-x — <Zi = p 3 , and n is the space-time dimension. 
The usual "causal" prescription for singularities in pseudo-Euclidean momentum 
space is understood, 

((q + k) 2 )- v ~{{q + kY+iQ)-», (2) 

When the powers of denominators and n are related by V\ + v-i + v 3 = n, a very 
simple result can be obtained for such a "unique" triangle |Q, ^] : 

J(n; v xM ) =W 2 » 1 -"n r(n/2 7' ) (^r-"/ 2 . (3) 
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If ^1 + 1/2 + v 3 = n — 1, then the following identity || (see also || [l(|) holds for 
such "semi- unique" triangles: 



{ ViJ(n; v\ + 1, v 2 ,v 3 ) +f2J(n; v u v 2 + 1, v 3 ) + v 3 J{n; v u z/ 2 , v 3 + 1)} | s „. =n _ 
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We shall use so-called "uniqueness" conditions (||) and (||) below, when evaluating 
the two-loop planar diagram. 

In applications the case v\ = Vi = v 3 = 1, n = 4 is important. Let us denote 

C^(pIpIpI) = J(4; 1,1,1). (5) 

Note that (||) can be represented in terms of a two-fold Mellin-Barnes integral fll| 

^ ioo ioo 

^ l ^ T —^—- [ [ dudvx u y v T 2 {-u)T 2 (-v)T 2 (l + u + v), (6) 
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where 

* = 4 > y = 4 , (7) 

and (here and below) the integration contours are chosen so as to separate the 
"right" and "left" series of poles of gamma functions in the integrand (see, e.g., 



[[12) ) . Analogous Mellin-Barnes representation for the case of arbitrary values of n 
and Vi can be found in Jl3], Oj. 

Using Feynman parameteric representation (or (0)) yields 
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with 
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The integral fj) can be evaluated in terms of dilogarithms Li 2 (z) (see, e.g., p5[), 



(a, y) = 2 (Li 2 (-px) + Li 2 (-py)) + \n(px) ln(py) 
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with 



X(x,y) = y/(l-x-y) 2 ~4xy 
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Note that representations of such type are well known for triangle diagrams (see, 
e.g., |l6[]). By simple dilogarithm transformations, formula (fw|) can be turned into 
the result obtained in [O . In the paper jPiJ some problems ofanalytic continuation 
of the function $W(x, y) ( |io| ) in the region of positive x and y were examined. If 
we consider, for example, negative values of both variables x and y, we should take 
into account the prescription (||) . This requires the following substitutions in ([H]) : 

ln(px) — » ln(— px) + iira ; hi(py) — > ln(— py) + z7T(j , (12) 



where a = sgn (p§) . 

Let us consider now a four-point "box" diagram in Fig. 2. All external momenta 



Figure 2: 

are taken to be ingoing (ki + fc 2 + k$ + ki = 0). We shall denote this diagram as 
D^[kl,kl,kl,kl,s,t) with 



s = (fci + k 2 ) 2 , t = (fc 2 + fc 3 ) 2 - 



(13) 



Then, by use of Feynman parametrization, uniqueness condition (^) and Mellin- 
Barnes representation, it is possible to show pS| that 
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where 



X = Mi , y^^i. (15) 
si si v ; 



Comparison with the representation (^) gives 

D ^ (k 1 , k 2 , k 3 , fc 4 , s, i) = C (fcj^^, fc 2 fc 4 , si), (16) 

or 

D«(^fci4fciM) = ^* (1 ^,n (17) 

s r 

with the same function $W as for triangle diagram. Thus, the result obtained 
for the box diagram contains only two dilogarithms. For negative x and y (this 
corresponds, e.g., to a physically interesting case when s and kf are positive while 
i is negative) , it is necessary to use a prescription of the type of ( |l2| ) (for details see 
injTTJ]). We note a useful fact of "pairing" of variables in the four-point function 
([16 fallowing reduction to the three-point function. Below we shall see that an 
analogous property also occurs for two-loop ladder diagrams. 

3. In this section we shall consider the main steps of evaluating two-loop thrce- 
and four-point ladder graphs. A more detailed information can be found in |fjq| . 

Let us consider first a three-point two-loop ladder (planar) diagram shown in 
Fig. 3 (pi + P2 + P3 = 0). The Feynman integral corresponding to this diagram 




Figure 3: 

(with unit powers of denominators) can be written as 

CW(pIpIpI) = J r2 ^ - f )2 \ p2 _ r)2 CU(( P1 + rf, ( P2 - r)\pl), (18) 

where is the one-loop function (||). 

To calculate , it is convenient to use the "uniqueness" method, by analogy 
with the paper JTjj (where propagator-type ladder diagrams have been examined). 
To do this, let us consider a special analytic regularization of this diagram (see 
Fig. 3), where we replace unit powers of denominators by (1 + Si) , provided that 
Si + 62 + S3 = 0. Applying relations (||) and (ji|) to this regularized diagram gives 
the following result (at n = 4): 
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where one-loop integrals J arc defined in ([j]). Note that integrals on the r.h.s. of 
(|19|) can be transformed by use of the formulae 



J(4; 1, f , 1 + S) = (pi)- d J(A; 1 + 8, 1 - 5, 1) = {p z 2 )- d J(4; 1 - 8, 1 + £, 1). (20) 

These relations can also be obtained from (|^) and (f|). 

By use of ( po|) and Feynman parameters, after some transformations one can 
obtain the following representation: 



.7(4; 1,1,1 + *) = 



(pi) 1+s s 
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(21) 



Inserting (21) into (E2h we get (as Si — > 0, J2$i — 0) 



where 
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(23) 



This integral can be easily calculated in terms of polylogarithms Li^ (z) (see [[L5|), 

l 



Li at (z) = 



(-1) 
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ln^C 



(24) 



So, we arrive at the following result for the two-loop three-point diagram of Fig. 3: 



( 2 ) (as, y) = i {e (Li 4 (-/«) + Li 4 (-py)) + 3 In | (Li 3 (-px) - Li 3 

+ \ ln2 1 ^ Li2 + Li2 + 1 1v2 (p x ) ln2 (py) 
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(25) 



where \{x, y) and p(x,y) are defined in (|l l|) . For negative values of x and y, one 
has to use (J12|) . 

Let us consider now the two- loop four-point diagram ( "double box" ) presented 
in Fig. 4. All notations correspond to the one-loop case (in particular, s and t are 
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Figure 4: 



denned by (|13|)). The corresponding integral can be represented as 
Z)' [k\ , fc 2 , fc 3 , fc 4 , s, t) = 



r 2 (/c3 + r) 2 (fc4 — r) 2 
x£>«(fc 2 , A; 2 , (fca + r) 2 , (fc 4 - r) 2 , s, (fc 2 + fc 3 + r) 2 ), (26) 

where is one-loop function (see Fig. 2). Note that the s-channel of the diagram 
in Fig. 4 corresponds to a "horizontal double box" with initial particles momenta k\ 
and &2' On the other hand, the t-channel corresponds to a "vertical double box" 
with initial momenta fc 2 and fc 3 . 

By using the same technique as for one-loop functions, one can construct four- 
fold Mellin-Barnes representations for and (see Comparing these 
representations we find 

D^\klkl,kl,kl,s,t)=tC^{klkl kjkl st), (27) 

or 

D^{klklktkts,t) = ^ ^(X,Y) (28) 

with X and Y defined by (15). We see that in two- loop case we also have obtained 
a "pairing" of four-point function arg uments. As a result, D 1 - 2 ^ can be reduced 
to three-point function (p^) -(|25T) (where we should make the substitutions 

p| — ► st, x — > X, y — » y ; sec (|15|)). Thus, formula ( p5|) (combined with ( p8| ) 
and (j^)) yields a representation of the "double box" i> 2 ) (Fig. 4) in terms of 
polylogarithms (Q). 

4. In this section we shall examine an interesting example of non-planar graphs: 
two-loop three-point crossed diagram [p\,p\,p^) shown in Fig. 5. We see that 
each single loop of this diagram corresponds to four-point function 

To evaluate the diagram in Fig. 5, it is convenient to use Fourier transform to 
coordinate space. It is easy to show that the "topology" of the diagram remains the 
same in the x-space (see Fig. 5): there are three points A±, A 2 , A3 (associated with 
external vertices), and there are two points V\ and V% connected by lines with each 
of Ai . Note that in the x-space both integrations (with respect to the positions of 
V\ and V2) are independent. Therefore, the crossed diagram in Fig. 5 factorizes 
and, using the x-space "uniqueness" condition (see, e.g., in ||) and returning to 
the momentum space, we find a simple relation: 

& 2 Hp\,pl,pl) = (c^Xp\,pI,pI)) 2 , (29) 




Figure 5: 



Figure 6: 



or (see (||)) 

C^iplplpl) = (^) 2 {^\x, V )) 2 (30) 

with Q^ix^y) defined by (p|)-(^0|). So, the result for the crossed diagram in Fig. 5 
is expressed in terms of the proHuct of dilogarithms. 

We can also use the formulae (p9|)-(|3"o|) to check the known result for three- loop 
crossed propagator-type diagram (see Fig. 6), 

B^\k 2 ) = j '^^(^(Hr)',^). (31) 



Representing in terms of a product of one- loop functions (P9|)-(p0[) and 
using parametric representation for <&W functions (this yields integrations over 
the parameters £ and £'), we obtain that the integral over r corresponds to a box 
diagram (Fig. 2) with ingoing momenta fci = (1 — &2 = — £(1 — £) _1 ^> 

fc 3 = £'(1 - £')- 1 fc, fc 4 = -(1 - £')~ lfc - The logarithms ln((fc + r) 2 /r 2 ) occurring 
in the numerator can be transformed into derivatives with respect to the powers of 
denominators by 



1 , i f(k + r) 2 \ ( d j 1 
■ In 



(k + r) 2 \ r 2 J \dSi ( r 2 ) 1+S ((k 



(32) 
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By use of the same technique as in Section 2 the box integral with shifted powers of 
denominators (1 + S, 1—5) can be reduced to the three-point function (J20|) . Using 
the parametric representation ( J2l| ) and making evident substitutions oFvariables, 
we arrive at the well-known result (see, e.g., pQ]): 



It should be noted that by using the representation ( |2l| ) the same result (^3|) can 
be obtained also for the planar three-loop propagator-type diagram. 

5. In the present paper we considered (via two-loop examples) the evaluation of 
planar and non-planar diagrams with three and four external lines with arbitrary 
momenta. We used Feynman parametrization, "uniqueness" conditions @-([|), 
Mellin-Barnes contour integrals and Fourier transformation. For the ladder dia- 
grams, it is shown that the corresponding four-point functions can be reduced to 
three-point ones (see (nw and (pTT)). Note that analogous formulae also occur for 
ladder diagrams with arbitrary number of rungs. 



The results (10) and ( [25| ) are presented in terms of polylogarithms with simple 
arguments (it is convenient to use such expressions in realistic calculations). Note 
that applying relations (||) and (|j) makes it possible to consider ladder diagrams 
with any numbers of loops. These results also can be expressed in terms of poly- 
logarithms (p4h . In the general case of L-loop ladder diagram, the highest order of 
occurring polylogarithms is equal to 2L. 

It should be noted that, if external momenta vanish, we get infrared (on-shell) 
singularities. For example, in the papers [^l[ [2^, |23| the diagrams of Fig. 3 and 
Fig. 5 have been considered at p\ = p\ = by use of dimensional regularization, and 
singularities have appeared as the poles in e — (4 — n)/2. The leading singularities 
of and 

are of the order of 1/e 4 (note that the coefficient at 1/e 4 in 
is four times as large than one in C' 2 '). In our (four-dimensional) approach in this 
case we can put p\ 2 = M 2 ( t 1 ~ * 0), and the singularities will be manifested as 
powers of In fx. It is easy to see from (|30|), ( |To| ) and (p2|), ( p5| ) that in this case the 
leading singularities of C' 2 ' and C' 2 ' are given by In 2 (pa;) ln 2 (py) ~ In /x, the ratio 
of coefficients being also equal to four. On the other hand, there is a problem how 
to relate more soft logarithmic and 1/e singularities. 
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